Classification of KPZQ and BDP models by 

■ multiaffine analysis 

o : 

^ ! HiROAKI KATSURAGI*and Haruo Honjo 

• Department of Applied Science for Electronics and Materials, 

C/^ ■ Interdisciplinary Graduate School of Engineering Sciences, 

'sj" ! Kyushu University, 6-1 Kasugakoen, Kasuga, Fukuoka 816-8580, Japan 

February 8, 2008 

Oh: 

d ■ 

• 1— I ' 

' Abstract 

' ' ' We argue difFerences between the Kardar-Parisi-Zhang with Quenched 

^ , disorder (KPZQ) and the Ballistic Deposition with Power-law noise (BDP) 

^ ' models, using the multiaffine analysis method. The KPZQ and the BDP 

I models show mono-affinity and multiaffinity, respectively. This difference 

. results from the different distribution types of neighbor-height differences 

' in growth paths. Exponential and power-law distributions are observed 

On , in the KPZQ and the BDP, respectively. In addition, we point out the 

. difference of profiles directly, i.e., although the surface profiles of both 

\l ' models and the growth path of the BDP model are rough, the growth 

, path of the KPZQ model is smooth. 

! 1 Introduction 

. ^ I The motion of rough surface has attracted much attention in recent decades 

■ ■ The most important feature is a dynamic scaling of surface width (root 
I mean square roughness) w. Consider a flat surface initially, and it grows to 

. - - 1 rough surface with time evolution. There are two scaling exponents a (roughness 

exponent) and /3 (growth exponent). The w is scaled by a scaling function ^' 
as w ~ x°''^{t/x^), where x, t, and z are the space length scale, time scale, 
and dynamic exponent z = a/ (3, respectively. The scaling function ^' consists 
of two parts depending on the argument u — t/x^\ "^{u) ^ for u ^ 1, and 
*(m) const, for u > 1 PE]. 

Since Kardar, Parisi, and Zhang (KPZ) proposed the continuous partial 
differential equation model '3', many theoretical and simulational studies have 
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Table 1: Experimental or observed results of a in anomalous scaling short range 
regime. All results show a > 1/2. 



Experiment 


a 


Reference 


Fluid flow 


0.73 


Rubio et al, 1989 ,4, 




0.81 


Horvath et al, 1991 [U 




0.65-0.91 


He et al, 1992 


Paper wetting 


0.63 


Buldyrev et al, 1992 


Bacteria growth 


0.78 


Vicsek et al, 1990 [HI 




0.78 


Wakita et al, 1997 "S* 


Burning front 


0.78 


Zhang et al, 1992 10^ 




0.81-0.89 


Myllys et al, 2000 11 


Crystal growth 


0.79-0.91 


Honjo and Ohta, 1994 12j 


Mountains profile 


0.73 


Matsushita and Ouchi,1989 [T3| 




0.68 


Katsuragi and Honjo, 1999 [T3| 



been done. The KPZ equation is written as follows, 

dth = vV'^h+^{Vhf +ri{x,t). (1) 

Where, h, v, X, and rj correspond to surface height, effective surface tension, lat- 
eral growth strength, and noise, respectively. They computed scaling exponents 
as a = 1/2 and (3=1/3 for 1-1-1 dimensional growing surfaces. Moreover, they 
found the scaling-law a + z = 2. This is a so-called KPZ universality scaling 
law. 

On the other hand, a lot of experiments of growing rough surface have been 
carried out. We show some typical results in Table As can be seen in Table 
n all experiments show a > 1/2 for anomalous scaling (short range) regime. 
This implies the KPZ equation model is not a sufficient model for growing rough 
surface phenomena. Short range anomalous scaling has been studied in other 
previous papers |15[ I16L [T7] . In ref. El rnany models for kinetic rough surface 
were considered. However, noise was assumed as Gaussian white in all cases. 
In refs. E| and El they used noises different from Gaussian white. We focus 
on this short range anomalous scaling regime from the point of view of the 
relation between multiaffinity and noise statistics, in this paper. Therefore, we 
performed numerical simulations of two variants of the KPZ model according 
to refs. 1161 andllTI In those models, the noise term is modified. In the one case, 
quenched disorder type noise is considered |16j. In the other case, uncorrelated 
power-law distributed noise is adopted . Both models can generate the large 
a{> 1/2) values. However, physical meanings of those two models are quite 
different. How can we distinguish those two models? The answer we present in 
this paper is the multiaffine analysis method. 

In this paper, we report on the result of numerical simulations of two variants 
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of the KPZ model. In the next section, multiaffinc analysis method is described. 
We introduce details of two models in fJ21 And we apply multiafhne analysis 
method to simulations results, in 21 In the roughness of growth path at each 
model is discussed. The comparison between our result and experimental data 
will be contained in this section. Furthermore, distribution forms of neighbor- 
height-differences are discussed, too. Finally, we conclude our results in the last 
section. 

2 Analysis methods 

The multiafhne analysis method is based on the multifractal concept ^ IllL I14L 
I18| . According to their results, q-th order height-height correlation function Cq 
is defined as, 

Cg{x,t) = {\Sh{x' ,t') - 6h{x' + x,t' + t)\'^)^,^t'- (2) 

Where Sh{x, t) is a deviation from mean height, i.e., Sh{x, t) = h{x, t) — {h{t))x- 
If we fix the temporal (or spacial) coordinate, we can define the q-th order 
roughness exponent aq and q-th order growth exponent (3q at the limit of a; 
and t — > as follows, 

Cg(x,0) (3a) 

CqiO^t)^^^^^. (3b) 

Obviously, aq and (3q with this definition quantify q-th order roughness scal- 
ing of surface profile 6h{x,t') (at fixed t') and growth path Sh{x',t) (at fixed 
x'), respectively. If the aq {/3q) varies with q, the surface profile (growth path) 
is multiaffinc (multi-growth). If they show constant values independent of q, 
it means mono-affinity and mono-growth property. More details about multi- 
affinity can be found in refs. ^ a-nd El The aq {(3q) just corresponds to 
the generalized dimension Dq in the multifractal formalism jl8| . When we find 
the changing Dq depending on q, we call that multifractal. Multiaffinity is an 
analogue of that. 

We can assume another definition of the growing exponent. That is based 
on the q-th order surface width Wq defined as follows, 

i=l 

where 6hi, L are height deviation at i-th site {Shi = hi — {hi)) and the number 
of sites, respectively. Then, the scaling forms can be written as Wq ~ a;"' and 
Wq ~ t^'' . We use this definition of (iq later. Using those Cq and Wq analysis 
methods, we classify two models of growing rough surface and discuss its original 
mechanism. 
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3 Models 



3.1 KPZQ model 

We examine two variants of KPZ model. The one is the Kardar-Parisi-Zhang 
with Quenched disorder (KPZQ) model that employs quenched disorder instead 
of time dependent noise. This model can be written as follows by a transform 
of KPZ noise term ri{x, t) 77(1:, h), 

dth = uV^h+^{Vhf + f + r]{x,h). (5) 

Where / is a positive driving force term. If there is not any such positive 
driving force, growth of surface shall stop due to the effect of negative quenched 
noise. Csahok et ai, calculated the characteristic exponents as a = 0.75 and 
/3 = 0.6, by the dimensional analysis method ^B]- These exponents satisfy the 
KPZ universality scaling-law a + z = 2. The discretized version of eq. (0) with 
simple single step method is written as, 

h{x,t + At) = h(x, t) + At(^{h{x -l,t) - 2h{x, t) + h{x + 1, <)} 

+ '^{h{x + l,t)- h{x - 1, t)Y +V + r]{x, [h{x, t)])) . (6) 

Where [h{x, t)] denotes integer part of h{x, t). We use quenched disorder rj{x, h) 
with a uniform distributed noise [—1, 1]. The parameters A, /, and At are taken 
as 1.0, 0.05, and 0.01 in all simulations. These values are close to pinning 
threshold. There is a crossover between short range anomalous scaling and 
asymptotic KPZ scaling, in this model (201 ■ We show a surface profile example 
of the KPZQ model in Fig. [Ha). 

3.2 BDP model 

The other model is the Ballistic Deposition with Power-law noise (BDP) model. 
Consider a flat substrate. Then, the noise particles are dropped from above 
randomly. They stick growing aggregate when they reach its nearest neighbor 
and next-nearest neighbor. This growing rule corresponds to ultra discrete 
version of KPZ equation The rule is written as, 

h(x,t + l) = \aa.y\h(x — l,t) + rj{x — l,t),h{x,t) + ri(x,t), 

h{x + l,t) + -qix + l,t)]. (7) 

Zhang obtained large a{> 1/2) value using the uncorrelated power-law dis- 
tributed noise with this model |T7]. Noise amplitude distributes as Pir]) ~ -j^+i) 
for rj > 1 and -P(ry) = for otherwise. Barabasi et ai, found the multiaffinity 
in this model . Buldyrev et ai, calculated the exponents as a = and 

P = 2ft^d '^ith mean field idea [23j in the region /i < 2d + 3. Where d is a spatial 
dimension of substrate. In addition. We have revealed the asymptotic behavior 
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Figure 1: Typical surface profiles of (a) KPZQ model and (b) BDP model. The 
parameters are taken as system size L = 4096 and calculation step T ~ 10^ 
steps in both models. 



of aq and Pq for large q regime |19| . Some experiments show the power- law 
distributed effective noise with ^ ~ 3 ^2 ■ ^ surface profile example of 
this model at fi = 3.0 (d = 1) is shown in Fig. ^b). In this case, a and /? 
become 0.75 and 0.6, respectively. These values coincide with those obtained 
from the KPZQ model as mentioned in t l3.1l Therefore, we cannot distinguish 
those two models only by the analysis of a and /?. Thus we use multiafRne 
analysis henceforth. 

4 Simulations and results 

First, we compute Cq of the KPZQ model. We show log-log plot of q-th root 
of Cq as a function of x (Fig. 12) and t (Fig.O. Each Figure displays 10 curves. 
Each curve corresponds to <? = 1, 2, . . . , 10, increasing from the bottom to the 
top. Fully parallel Cq curves are in small x (or t) regime. This implies the mono- 
affinity of the KPZQ model. The slope of curves in Fig. |21 indicates aq = 0.75 
in short range scale. This value agrees with the theoretically predicted value 
a = 0.75. We use the system size L — 4096, calculation steps T — 10^ steps, 
and ensemble number N — 100. 

Similarly, curves in Fig. |21 have the same slope. It means mono-growth 
property of the KPZQ growth. However, its value in the short range scale 
Pq = 1.0 is not trivial. While Csahok et ai, calculated the KPZ universality 
compatible result (3 — 0.6, we obtain (3 — 1.0. We use L = 1024, T = 10'' steps, 
and iV = 1 for this calculation. In Figs. |2] and |21 the KPZ scaling {aq = 0.5 and 
Pq — 1/3) can be slightly observed in the long range scale. 

We can also evaluate Pq using the Wq scaling according to eq. Q . We show 
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Figure 2: |logCq(x,t') as a function of logx of the KPZQ model. Each 
curve has the same slope 0.75 that indicates mono-afHnity of the KPZQ model. 
Moreover, it is consistent with the dimensional analysis of Csahok et al. The 
parameters are taken as L = 4096, T = 10^ steps, and N = 100. 
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Figure 3: ^ log (a;' , as a function of logt of the KPZQ model. The slope 
1.0 conflicts with the dimensional analysis of Csahok et al.. This means that the 
KPZ universality is destroyed. The parameters are taken as -L = 1024, T = 10^ 
steps, and N = 1. 



-2-1012 
logio ' 

Figure 4: logit;^ vs. logt of the KPZQ model. Although there include a little 
fluctuations, we can confirm the fit to the slope 0.6. The parameters are taken 
as L = 1024, T 5 X lO"*, and N = 200. 



the log-log plot of Wq as a function of t in Fig. ^ The curves correspond to 
q = 1, 2, . . . , 10 from the top to the bottom. The slope corresponds to Pq{— 0.6). 
Mono-growth property can be seen again, and the value agrees with theoretical 
result (the KPZ universality is recovered). 

In addition, we introduce one more quantity in order to describe the growing 
rough surfaces. We consider a growth path dh{x',t) at fixed x' , and define a 
neighbor-height-difference (or local growth velocity) rjt as follows, 

r]t = \5h{x',t)~6h{x,t+l)\. (8) 

The size distribution of r/t is displayed in Fig. [3 As can be seen in Fig. El the 
distribution has the exponential tail. 

We have reported the results on the BDP model in our previous paper jl9| . 
Let me summarize the result here. The BDP model exhibits multiaffinity and 
multi-growth property, i.e., aq and f3q varies with q. The size distribution of 
neighbor-height-differences has the power-law tail. This power-law allows us to 
use the multifractal dimensional analysis method in order to obtain the asymp- 
totic function of Uq and (3q for large q regime {1/q < aq < 2/{q + 1) and 
f3q ~ 1/q). These behavior of aq and (3q are quite different from those of the 
KPZQ model. As mentioned before, the a, and (3q become constants in the 
KPZQ model. Note that the size distribution of neighbor-height-differences in 
the KPZQ model has the exponential tail. Therefore, we cannot apply the 
multifractal dimensional analysis method to the KPZQ model. We have also 
reported on the mono-affinity of the KPZQ model in ref. 1251 however, it was 
not performed in sufficient system size. The results could not be quantitatively 
reliable enough. 

As described above, the KPZQ model shows mono-affinity and the BDP 
model shows multiaffinity, indeed. And those relate to the distribution types 
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Figure 5: Distribution of neighbor-height-differences of the KPZQ model. Clear 
exponential part about P{r]t) ~ exp(— lOyyt) can be observed. The parameters 
are taken as L ^ 1024, T = 10"* steps, and N = 10. 



of neighbor- height-differences; exponential or power-law, respectively. That is, 
we can distinguish the very similar profiles, {e.g., Fig. ^ using the multiafRne 
analysis. 

5 Discussion 

In the KPZQ model, q-th order growth exponent /3g depends on measuring 
methods. The Cq{x',t) based method yields Pq = 1.0. Contrary, the Wq method 
provides Pq = 0.6 that concurs with theory and satisfies KPZ universality. Be- 
cause Cq{x',t) expresses the self-affinity of growth path, the box-count fractal 
dimension Db oi the growth path can be written as Z?s = 2 — /3i |14[I18| . In the 
present case, we obtain jSq = 1.0 for all q, therefore, Db becomes 1.0. This infer 
that the growth path is not a rough curve, but a smooth 1-dimensional one. We 
show examples of growth path of each model in Fig. As expected, smooth 
and rough growth paths are confirmed for the KPZQ and the BDP model, re- 
spectively. In this mean, the value of Pq of the BDP model does not depend 
on measuring methods. The growth paths of the BDP model are ordinal rough 
profiles [3 El- Moreover, the scale of 6h in Fig. (S^a) is significantly smaller 
than that of Fig.El^b). These data indicate that the growth of the KPZQ model 
is calmer than that of the BDP model, at least in the small t regime. 

Meakin pointed out the Hurst exponent and the roughness exponent are not 
always same [3^. He has discriminated self- affinity and roughness scaling. His 
thought was limited in spatial dimension. Our result on Pq corresponds to the 
discrimination between the self-affinity of growth path and the growth scaling. 
The Cq method provides the former, and Wq mehod provides the latter. More- 
over, slopes for long range scale in Figs.|21and|31seem to be calm. This probably 
implies the asymptotic KPZ behavior (a = 0.5 and /3 = 1/3). Simulations in 
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Figure 6: Typical growth path of (a) KPZQ model and (b) BDP model. The 
parameters are taken as L = 1024 and T — 10^ steps for both models, and 
H = 3.0 for the BDP model. 



larger system are needed for more detail analyses. 

Slow combustion front of paper shows multiafhnity jllj . The size distribution 
of neighbor-hcight-differences (effective noise) of this system obeys power-law. 
These characteristics are consistent with the BDP model results. On the other 
hand, the profile of the Hida mountains shows mono-affinity ^3]. Furthermore, 
the size distribution of neighbor-height-differences of this profile shows almost 
exponential form. These features follow the KPZQ model results. However, 
rough surface of paper combustion seems to be caused by quenched disorder 
on papers. And mountain profiles probably result from the rain-fall erosion. 
Thus, physical mechanism of each model seems to pass each other, intuitively. 
Recently, Maunuksela et al., and Myllys et al., performed numerical simulations 
based on real paper noise |27[ I28| . They computed the value of coefficients in 
the KPZ equation from experiments, and used the quenched noise measured 
from the real paper. They obtained mono-scaling result for this type of realistic 
quenched disorder simulations. We could not identify the specific origin of each 
noise type perfectly. More comprehensive analysis about experimental data is 
necessary to understand universality class in detail. This problem is an open 
question. There are many other kinetic rough surface models 1 . It must be 
interesting to study on these models from the viewpoint of multiaffinity. 

We have proved the validity of multiaffine analysis method in this paper. 
We believe that this analysis plays an important role to classify other scaling 
phenomena. A weak point of multiaffine analysis method is that it is sensitive 
to fluctuations of data compared with single roughness and growth exponent 
analysis, particularly to calculate higher order moment (large q regime). While 
the method is powerful, it requires numerous data. The improvement of the 
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method and application to other physical systems are future problems. 



6 Concluding remarks 

Wc applied multiafRne analysis method to the KPZQ and the BDP model, 
and found some features as following. The KPZQ and the BDP models have 
multiaffinity and mono-affinity, respectively. The size distribution of neighbor- 
hcight-diffcrcnccs obeys exponential form in the KPZQ model, and obeys power- 
law form in the BDP model. While growth path of the BDP model is rough, 
that of the KPZQ model is smooth. Although the KPZQ and the BDP models 
cannot be distinguished only by measuring a and /3, multiaffine analysis makes 
these differences clear. 
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